First order phase transitions in Ferromagnet/Superconductor layered structures 
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We study the thermodynamics of clean structures composed of superconductor (S) and ferromag- 
net (F) layers and consisting of one or more SFS junctions. We use fully self consistent numerical 
methods to compute the condensation free energies of the possible order parameter configurations 
as a function of temperature T. As T varies, we find that there are phase transitions between states 
characterized by different junction configurations (denoted as "0" or 'V according to the phase 
difference of the order parameter in consecutive S layers). We show that these transitions are of 
first order. We calculate the associated latent heats and find them to be measurable. 
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A plethora of new ideas and devices has been emerg- 
ing from the study of nanostructures as they pertain to 
the field of spintronics 1 . An important part of this devel- 
opment has occurred through the study of the rich and 
varied phenomena that occurf, in heterostructures involv- 
ing superconductors (S) and ferromagnets (F). 

The physics of such F/S heterostructures is dominated 
by the proximity effects that arise from the competition 
between superconducting and magnetic orderings in the 
materials comprising the structure, with each of the cor- 
responding order parameters penetrating into the other 
material. These effects follow from normal and Andreev* 
reflection processes at the interfaces. In the latter pro- 
cess, an electron encounters one of the interfaces, is con- 
verted into a hole with opposite spin and traverses in 
the opposite direction. For S/N interfaces (where N is a 
non-magnetic, non-superconducting metal) the dynamics 
of charge transport is degenerate with respect to the elec- 
tron spin quantum variable. This is not true, however, in 
the case we consider here, where superconducting regions 
are separated by magnetic interlayers. 

In SFS trilayers, as well as in multilayers built from 
a sequence of such structures (SFSFS. ..), the spin- 
splitting effect of the magnet produces important and 
nontrivial changes in Andreev and other scattering pro- 
cesses. The exchange field in the ferromagnet breaks the 
time-reversal symmetry and generates a superconducting 
state where the Cooper pairs acquire a finite momentum 
resulting in a spatial modulation of the superconducting 
pair amplitude in the F region*. Depending on the geo- 
metric and material characteristics of the SFS trilaycr, its 
thermodynamic equilibrium state can be a "0" or a 'V 
state, depending on the value of the phase difference be- 
tween the superconducting order parameters in the two 
S electrodes. It is this twofold possibility that lies at the 
foundation of the many spin-based switching phenom- 
ena, which in turn are the basis for devices, including 
superconducting 7r qubits* and memory elements*. For 
larger SFSFS. . . S type heterostructures the order param- 



eter may or may not flip between any pair of consecutive 
S layers, leading to a variety of possible configurations, 
which can be characterized as a sequence of and ir junc- 
tions. 

Continual advances in nanoprocessing methods have 
made it possible to fabricate high quality structures con- 
taining SFS junctions, which have encouraged further 
study of these systems. From the thermodynamics point 
of view, a to 7r transition as the temperature was varied 
was inferred**^** from Josephson current measurements in 
Nb/CuNi/Nb junctions. The Josephson coupling in simi- 
lar structures was also found to cross over from positive to 
negative, depending on the F layer thickness^, indicating 
again a to 7r switch. Under many experimental condi- 
tions, a change in the second-order Josephson coupling 
component was associated^*** with the — it crossover 
while for other conditions* the nonlinear current-phase 
relation did not reveal any change in the second-order 
Josephson coupling. 

The thermodynamics of F/S structures has been di- 
rectly examined largely in terms of studying the tran- 
sition temperature to the normal state*i*i*******. A 
first order phase transition to the normal state has been 
predicted^* in spin-valve ferromagnet-superconductor- 
ferromagnet (FSF) nanostructures with parallel magne- 
tization alignments in thin F layers. For antiparallcl ori- 
entation of the magnetization, the transition was found 
to be always second order—. The application of a bias 
voltage 13 up to a critical value can cause the phase tran- 
sition to be first order for both P and AP alignments 
of the magnetization. For F/S bilayers with stripe do- 
main structure*^, superconductivity appears via a first 
order transition. In some cases, with the action of an 
applied field, the compensating screening currents can 
enhanesi* the superconducting state. The theoretical 
literature that examines the thermodynamic transitions 
between and ir configurations in SFS type structures 
is more sparse. It has been arguecU* that the transi- 
tion that takes place in Josephson junctions is continu- 
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ous when the F layer thickness is uniform but it rounds 
off when it is variable. To further understand the un- 
derlying competition between the various possible states 
and to better tailor these structures for practical appli- 
cations, it is imperative to clarify the thermodynamics of 
systems potentially containing 7r junctions by investigat- 
ing the parameters that may influence a first or second 
order phase transition from a state to a 7r state or vice 
versa. 

The objective of this paper is to clarify some of these 
issues by rigorously considering the thermodynamics of 
clean layered systems consisting of one or more SFS junc- 
tions so as to identify and characterize any phase tran- 
sitions involving — tt flipping. It has been shown at 
low-temperatures^ that for given S and F widths, ex- 
change energy and other material parameters, multiple 
spatial configurations of the self consistent pair ampli- 
tude can exist as local minima of the energy. The larger 
the number of layers, the more combinations were found 
to be possible. Among the various solutions, the ground 
state was found from accurate condensation energy com- 
putations. Here on the other hand, we have the more 
ambitious objective of studying the possible competing 
states as a function of temperature T through a careful 
analysis of the free energy differences. We find that, as T 
varies, phase transitions associated with flipping of SFS 
junctions from a to a 7r state occur and that there is 
a discontinuity in the entropy at such transitions, which 
therefore are of first order. We calculate the correspond- 
ing latent heat, and find that its magnitude is observable 
given current experimental capabilities. 

We study layered S/F systems containing a number 
Nj of SFS junctions. We consider in particular Nj = 1, 
the important case of a single junction, and, to show the 
richness and variety of the possible outcomes, also the 
case Nj = 3. We compare as a function of T the con- 
densation free energies AJF(T), defined as the difference 
AJF(T) = T S {T) - T N {T) between the free energies of 
the normal and superconductor states of each of the sev- 
eral possible self-consistent configurations. It is in prin- 
ciple very difficult to compute condensation free energies 
to the required accuracy: one needs only to recall that 
in the bulk and at zero T, AT equals^ -(l/2)JV(0)Ag 
(A is the bulk gap). Since Fn is much larger, of order 
N(0)uip (ujd is the Debye frequency), AT is the small 
difference between two much larger quantities. Thus, one 
is faced with the tough numerical task of computing two 
quantities accurately enough so that the small difference 
between them can be reliably established. The problem 
is exacerbated for our systems because, as will be seen, 
the difference between the condensation free energies of 
the various possible states (e.g., the and tt cases for 
SFS trilayers) is at best a small fraction of their individ- 
ual condensation free energies. The numerical methods 
we use do resolve these small differences. These calcula- 
tions are carried out using the spectra and pair potentials 
obtained by fully self consistent numerical solution of the 
Bogoliubov deGennesi^ (BdG) equations. 



The systems studied here consist then of either 3 or 7 
layers, containing 1 or 3 SFS junctions respectively. We 
denote the thickness of the F layers by <ip and that of 
the S layers by ds and assume the interfaces are sharp. 
The procedures that we use to numerically solve the BdG 
equations for a clean syst em in a fully self consistent man- 
ner are detailed in Refs. Il7l20l and the details need not 
be repeated here. As explained there, for many relevant 
values of the geometric and material parameters several 
self consistent solutions, that is, local minima of the free 
energy, often coexist at T — 0. For example, solutions 17 
of both the and tt type may be possible for one junc- 
tion, while for the three junction case more than one of 
the possible symmetric states (000, iririr, irOir and 07r0 
in the obvious notation denoting the state of each junc- 
tion) may yield a self-consistent solution to the problem. 
In such cases, the equilibrium state had to be found by 
comparing the respective energies. It was shown that the 
symmetry of the stable solution at T = can changed 
as one varies parameters such as dp, the interface barrier 
height, the exchange field of the magnet, or the Fermi 
wavevector mismatch between F and S materials. The 
switchover from one stable state to another occurs as the 
corresponding energies cross at certain points in parame- 
ter space: "phase transitions" take place as a function of 
such parameters. They may occur also in the dirty limit 1 '' 
as a function of dp. The issue is how these phenomena 
are related to transitions found experimentallji2i2ii£ at 
finite T. 

Here we investigate the true thermodynamics, that is, 
how the equilibrium state of a given system depends on 
T, and the nature of the corresponding phase transi- 
tions. When self-consistent numerical solutions of the 
BdG equations are possible for multiple order parameter 
configurations, we evaluate the condensation free ener- 
gies of the competing states. For this purpose we use the 
method of Ref. l2ll which requires only the energy spec- 
trum and the position dependent order parameter, but 
not the eigenfunctions. We find that in many cases there 
are first order phase transitions as a function of T, and 
that the latent heats are quite appreciable. 

To locate phase transitions as a function of T, we 
searched near regions of parameter space where crossings 
occur— at zero T. We focussed on the Fermi wavevector 
mismatch parameter defined as A = Efm / Eps, where 20 
Eps is the Fermi energy (bandwidth) in the S material, 
while in F we have E± = Efm ± h for the majority (+) 
and minority (— ) bands. Here h is the exchange field. 
We also considered the geometric parameter, dp . For 
Nj = 1 we take ds equal to the zero temperature coher- 
ence length, £o , while for Nj = 3 we double the thickness 
of the two inner S layers to 2£o- This mimics a more 
symmetric situation, since the outer layers are in contact 
with only one F layer. We measure energies in units of 
EpSj lengths in units of kps> the Fermi wavevector for 
S, and set kps£,o = 100. The strength of the magnet 
is chosen as h = 0.2Efm throughout. We assume that 
all magnetic layers are aligned and that there is no ox- 
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FIG. 1: (Color online) Normalized condensation free energy 
(see text) as a function of reduced temperature. The 
results in the top panel are for a one junction (SFS) system 
at A = 0.535 and ^fs^f = 7.5. Results for the two possi- 
ble order parameter configurations are plotted as indicated. 
Those in the bottom panel are for a three junction system 
at A = 0.45 and ^fs^f = 10. Only the two lowest compet- 
ing free energy configurations are plotted. The insets display 
the normalized pair amplitude F(Z) at the transition points, 
where Z = &fsz. 

ide inter layer barrier. We performed several checks of 
our numerical methods. For a system consisting of only 
one S layer with ds 3> £oj we quantitatively recover the 
textbook results^ for the thermodynamics, including the 
second order transition at the correct bulk value Tj? and 
associated specific heat discontinuity. 

Proceeding now to the results for Nj = 1, 3, the main 
plots in Fig. ^ show the condensation free energy AT, 
as defined above, normalized by N(Q)A%, which is twice 
the value of the condensation energy for a bulk S mate- 
rial at T = 0. In the top panel, AT is plotted versus 
reduced temperature T/T® for both the and ir configu- 
rations of a 1 junction system. These results correspond 
to kFsd-F = 7.5 and A = 0.535; similar results are found 
at this thickness for a range of A. Points at 0.01 intervals 
in the horizontal temperature scale are plotted, joined by 
straight segments. The slopes of the AT curves approach 
zero as T — ► 0, which indicates zero entropy at T = 0. 
The slopes can also be seen to approach zero as AT — ► 0: 
thus we find that the transition to the normal state is 
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FIG. 2: (Color online) Dimensionless (see text) condensation 
entropy S(T) (right scales) and energy U(T) (left scales). Top 
and bottom panels correspond to the same cases as in Fig. 
The vertical arrows locate the first-order transitions. 

second order, occurring at T c < T°. One can also see 
in the upper plot that the results for the and tt states 
cross at T x /T® as 0.28, with the tt configuration being the 
equilibrium state below T x and the state above. The 
difference in the slopes at T x represents a latent heat, dis- 
cussed below. Thus we predict that there is a transition 
and that it is first order. The inset displays the spatial 
profile of the corresponding pair amplitudes at T x . The 
quantity plotted, F(Z), is the pair amplitude as usually 
defined^, normalized to its bulk S value at T = 0. Z 
is the dimensionless distance Z = kp$z normal to the 
layers. A discontinuity in the absolute value of F(Z) at 
the transition can be noted. 

The lower panel of Fig. ^ shows AT(T) similarly plot- 
ted for a 3 junction (7 layer) system with IcfscIf = 10 
and A = 0.45. Points are again taken at intervals of 
0.01 on the horizontal scale. In this case only two (7T7r7r 
and 7r07r) of the possible configuration states compete^ 
as candidates for lowest free energy over the range of T 
studied. The other possible states are therefore irrele- 
vant and omitted. The curves in the lower panel display 
the same characteristics as in the 1 junction case: zero 
first derivatives as T — > 0, a second order transition to 
the normal state, and a first order transition at an inter- 
mediate T. In this case the transition is at T x /T® 0.27 
and from a 7r07r state at lower T to tttttt at higher T. 
The inset again shows the corresponding F(Z) at the 
transition. Thus in this more complicated example, the 
transition involves the inner junction flipping to tt as 
T is increased. Again, such first-order transitions exist 
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FIG. 3: (Color online) Variation of the latent heat with T x . 
The entropy change, L, at T x , normalized to the normal state 
specific heat at T®, is plotted as a function of T x /T° for the 
one and three junction cases. 



over a range of A at fixed dp. 

We can now derive the entire thermodynamics. Thus, 
Fig. [21 shows the dimensionless condensation entropies, 
S(T), obtained by differentiating the results for AT(T) 
(Fig. ^| with respect to T/T®. We show also the corre- 
sponding condensation energies, U{T), computed from 
standard thermodynamic relations, normalized in the 
same way as AT. In calculating the derivatives, a poly- 
nomial form was fit to the AT data in Fig.^ The values 
found were equivalent to those obtained by taking a dis- 
crete derivative of the AT curves and fitting a polyno- 
mial to the resultant data. The top and bottom panels 
of Fig.|21correspond, respectively, to the 1 and 3 junction 
cases in Fig. ^ As mentioned previously, the entropies 
go to zero smoothly as T — > and T — > T c . On the 
other hand, it is quite obvious that the transitions at T x 
(vertical arrows) are indeed first order: the entropies of 
the two states involved are not equal at that tempera- 
ture. The condensation energy curves resemble those for 
a bulk superconductor. By taking a further derivative, 
the specific heat can also be obtained. For any given 
state the quantities S(T), U(T) and AT all go to zero 
at the same temperature, which is the computed value 
of T c . The energies and entropies cross at temperatures 
above T x : one can see that both entropy and energy play 
important and subtle roles in the first-order transition. 



As one varies A, T x moves up or down smoothly and 
monotonically until the transition disappears when T x 
reaches either or T c . The difference between the en- 
tropies at T — T x , which measures the latent heat, L, is 
quite appreciable and does not depend drastically on T x 
(nor equivalcntly on A) over most of the range. These 
points are illustrated in Fig. [31 where L is shown for 
both the 1 and 3 junction systems. We show L (calcu- 
lated as discussed in connection with Fig. 0) at several 
T x , obtained by varying A. The sign of L is defined by 
subtracting the entropy of the stable state below T x from 
that above T x . So that the size of the effect can be ap- 
preciated, L is normalized to the specific heat at T* of 
a normal S material of the same volume. This is ap- 
propriate since for a normal metal the specific heat at 
T® equals its entropy. One can see that, for both one 
and three junction samples, the entropy jumps at T x can 
exceed 1% of the entropy present in a uniform bulk S 
sample at Tj?. 

Translating these results for the normalized measure of 
the latent heat into actual units for typical samples*, we 
find that they are of order of up to 10 w k B T or 1CU 13 J. 
Standard ac calorimetry techniques offer a resolution22i2^ 
at least one order of magnitude smaller. Indeed, even at- 
tojoule calorimetry has been achieved in electronic sys- 
tems, although by using multiple samples^ Our rigorous 
results are clearly larger than estimates madei^ for dirty 
systems. These estimates, based on partial consideration 
of the Josephson energy only, are in the nature of lower 
bounds. 

In summary, we have rigorously studied the thermody- 
namics of clean single and triple SFS junction systems. 
We have shown that as T is varied a given junction can 
flip from the state to the it state. The resulting phase 
transition is fi rst-o rder, in agreement with the experi- 
ments of Refs. I7l8lll0l for systems with sharp interfaces. 
The associated latent heats were found to be within cur- 
rent experimental resolution. The results presented here 
should be applicable over a broad range of material pa- 
rameters assuming ds does not exceed several £o- 
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